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Abstract. - We demonstrate possible scenarios for production of pure spin current and large tun- 
nelling magnetoresistance ratios from elastic co-tunnelling and crossed Andreev reflection across 
a superconducting junction comprising of normal metal-superconductor-normal metal, where, the 
normal metal is a one-dimensional interacting quantum wire. We show that there are fixed points 
in the theory which correspond to the case of pure spin current. We analyze the influence of 
electron-electron interaction and see how it stabilizes or de-stabilizes the production of pure spin 
current. These fixed points can be of direct experimental relevance for spintronics application 
of normal metal-superconductor-normal metal junctions of one- dimensional quantum wires. We 
also calculate the power law temperature dependence of the crossed Andreev reflection enhanced 
tunnelling magnetoresistance ratio for the normal metal-superconductor-normal metal junction. 
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. Introduction. — Two fundamental degrees of free- 
dom associated with an electron that are of direct interest 
to condensed matter physics are its charge and spin. Until 
^evy recently, all conventional clcctron-bascd devices have 
been solely based upon the utilization and manipulation 
of the charge degree of freedom of an electron. However, 
the realization of the fact that devices based on the spin 
degree of freedom can be almost dissipation-less and with 
very fast switching times, has led to an upsurge in research 
dctivity in this direction in recent years [1-3]. The first 
step towards realization of spin-based electronics (spin- 
tronics) would be to produce pure spin current (SC). From 
a purely theoretical point of view, it is straight forward to 
define a charge current as a product of local charge den- 
sity with the charge velocity, but such a definition cannot 
be straight forwardly extended to the case of SC. This is 
because both spin S and velocity v are vector quantities 
and hence the product of two such vectors will be a tensor. 

In this letter, we adopt the simple minded definition 
of SC, which is commonly used [4]. It is just the prod- 
uct of the local spin polarization density associated with 
the electron or hole, (a scalar s which is either positive 
for up-spin or negative for down-spin) and its velocity [4] . 
The most obvious scenario in which one can generate a 
pure SC in the sense defined above would be to have (a) 



an equal and opposite flow of identically spin-polarized 
electrons through a channel, such that the net charge 
current through the channel is nullified leaving behind a 
pure SC, or (h) alternatively, an equal flow of identically 
spin polarized electrons and holes in the same direction 
through a channel giving rise to pure SC with perfect can- 
cellation of charge current. In this letter, we explore the 
second possibility for generating pure SC using a normal 
metal— superconductor— normal metal (NSN) junction. 

Proposed device and its theoretical modelling. — 

The configuration we have in mind for the production of 
pure SC is shown in Fig. 1. The idea is to induce a pair 
potential in a small region of a quantum wire (QW) by de- 
positing a superconducting strip on top of the wire (which 
may be, for instance, a carbon nanotube) due to proximity 
effects. If the strip width on the wire is of the order of the 
phase coherence length of the superconductor, then both 
direct electron to electron co-tunnelling as well as crossed 
Andreev electron to hole tunnelling can occur across the 
superconducting region [5]. It is worth pointing out that 
in the case of a singlet superconductor, which is the case 
we consider, both the tunnelling processes will conserve 
spin. In order to describe the mode of operation of the 
device (see Fig. 1), we first assume that the S'-matrix rep- 
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Fig. 1: A 1-D quantum wire (carbon nanotubc) connected to a ferromagnetic (F) lead on the left and a normal (N) lead on the right. 
The thin strip in the middle of the wire depicts a 2-D layer of superconducting material deposited on top of the wire. 



resenting the NSN junction described above respects par- 
ity symmetry about the junction, particle-hole symmetry 
and spin-rotation symmetry. Considering all the symme- 
tries, we can describe the superconducting junction con- 
necting the two half wires by an S'-matrix with only four 
independent parameters namely, (i) the normal refection 
amplitude (r) for e (h), (ii) the transmission amplitude [t) 
for e (h), (Hi) the Andreev reflection (AR) amplitude (r^) 
for e (h), and (iv) the crossed Andreev reflection (CAR) 
amplitude {Ia) for e (h). If we inject spin polarized elec- 
tron (t e) from the left QW using a ferromagnetic contact 
and tune the junction parameters such that t and ta are 
equal to each other, it will lead to a pure SC flowing in 
the right QW (see Fig. 1). This is so because, on an av- 
erage, an equal number of electrons (t e) (direct electron 
to electron tunnelling) and holes (t h) (crossed Andreev 
electron to hole tunnelling) are injected from the left wire 
to the right wire resulting in production of pure SC in the 
right wire. Note that spin up holes (t h) implies a Fermi 
sea with an absence of spin down electron (which is what 
is needed for the incident electron (t e) to form a Cooper 
pair and enter the singlet superconductor). 

Renormalisation Group study. Wc now include 
the effects of inter-electron interactions on the S'-matrix 
using the renormalisation group (RG) method introduced 
in Ref . [6] . This was furthur generalized to the case of mul- 
tiple wires in Refs. [7,8] and to the case of 1-D NS junc- 
tion [9-11]. The basic idea of the method is as follows. The 
presence of scattering (reflection) induces Friedel oscilla- 
tions in the density of non-interacting electrons. Within 
a mean field picture for a weakly interacting electron gas. 
the electron not only scatters off the potential barrier but 
also scatters off these density oscillations with an ampli- 
tude proportional to the interaction strength. Hence by 
calculating the total reflection amplitude due to scattering 
from the scalar potential scatterer and from the Friedel 
oscillations created by the scatterer, we can include the 
effect of electron-electron interaction in calculating trans- 
port. This approach can be generalized to junctions of 



one-dimensional (1-D) QW with a single superconductor. 
In this case, there will be non-zero AR in the bulk of the 
wire due to proximity induced pair potential, besides the 
AR right at the NSN junction which turns an incoming 
electron into an outgoing hole. 

The RG equations for a NS case have been obtained ear- 
lier using bosonization [12-14] and using WIRG [9, 11]. In 
this letter we extend these WIRG results to the NSN case. 
To obtain the RG equations in the presence of AR and 
CAR for the NSN junction, we follow a procedure simi- 
lar to that used in Ref. [7]. The fcrmion fields expanded 
around the left and right Fermi momenta on each wire can 
be written as, V'is {x) = ,s (cc) e* ^ ^'o is {x) e"* ^ ; 
where i is the wire index, s is the spin index which can be 
t,i and I{0) stand for incoming (outgoing) fields. Note 
that ^'/(o)(x) arc slowly varying fields on the scale of fc^^^. 
For a momentum in the vicinity of fc^?, the incoming and 
outgoing fields can be Fourier expanded as: 



*fcs(a;) 



dk 



gj (fe + ^f) a; _|_ ^i{-k + kp) X 



(1) 



where 6^,^ is the electron destruction operator and d^^ is 
the hole destruction operator and we have allowed for non- 
conservation of charge due to the proximity effect. We 
allow for short-range density-density interactions between 

the fermions, Hint ~ 2 I dx dy p^gV{x — y) p^^, where the 
sum over the spin indices is assumed. 

Then the effective Hamiltonian, can be derived using a 
Hartree— Fock (HF) decomposition of the interaction. The 
charge conserving HF decomposition leads to the interac- 
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Fig. 2: The variation of |tp {=|iAp) is plotted as a function of 
the dimensionless parameter / where I = ln(L/d) and L is either 
Lx = hvp/kgT at zero bias or Ly = hvp/eV at zero temperature 
and d is the short distance cut-off for the RG fiow. The three curves 
correspond to three different values of V{0) and V{2kp) for the 
NSN junction. This case corresponds to the 5-matrix given by Si. 



tion Hamiltonian (normal) of the form 
-i{g2 - 2gi) f°° dx 



njN 

'lint — 



47r 



a 



(2) 



(We have assumed spin-symmetry and used r^^ = 7-^ = 
Ti.) This has been derived earlier in Ref. [7]. We use the 
same method, but now, we also allow for a charge non- 
conserving HF decomposition and arrive at the (Andreev) 
Hamiltonian 



ri\nt — 



-*(.9i + .92) 



dx 



(3) 



Note that, even if this appears non-charge conserving, 
charge conservation is taken care of by the 2e charge that 
flows into the superconductor due to proximity induced 
Cooper pairing. 

In Ref. [7], the perturbatively calculated correction to 
the reflection amplitude under exp [—iTi.^^ t] was derived 
to first order in a. For electrons with spin incident with 
momentum k with respect to kp, this was shown to be 



a short distance cutoff and 171 = V{2kF), 52 = ^(0)- It 
is worth noting that both gi and g2 scale under RG [6]. 
Hence, the values given above for gi and 172 are the micro- 
scopic values it takes in the original Hamiltonian defined 
at the short distance scale. 

Analogously, here, we calculate the amplitude to go 
from an incoming electron wave to an outgoing hole wave 

a' TAt 



under exp [—iTi^^ t 

{91 + g2)/27rfiu_F 



It is given by 



■ In(fcd) where 



Q^' ~ (91 + g2) T'^T^tWF ■ Note that the spin of the outgo- 
ing hole is always the same as the spin of the incoming 
electron, since the Andreev Hamiltonian also conserves 
spin for a singlet superconductor. We see that there is a 
logarithmic singularity at the fc ^ limit which implies 
that the lowest order perturbation theory is not enough to 
calculate correction to the reflection and AR amplitudes 
when the momenta of the incident particles are very close 
to the Fermi wave vector. Following Yue et al. [6], we sum 
up these most divergent processes using the "poor man's 
scaling" approach [15] to obtain RG equations for the nor- 
mal refection amplitude (r), the transmission amplitude 
(i), the AR amplitude (r^), and the CAR amplitude (f^) 
which are as follows, 



dr 

m 

drA 
dl 

dt 
Jl 

dtA 
dl 



r\+t\) -r{l-\r\^)) 



~a' {r|rAp + r\tAt] 

- a {\r\^rA + ttAr*] + y {ta 

- a {\r\^t + r*rAtA] 
-a' {\rA\h + rr\tA] 

- a {rV^i + \r\'^tA] 
-a {rtr\ + \rA\^tA] 



(4) 



(5) 



(6) 



(7) 



Note that when t — tA = ^ these equations reduce to the 
RG equations obtained in Ref. [9] for the case of NS junc- 
tion. 

Results and Discussion. — We propose two possible 
S*- matrices {Si and S2) that can be realized within our set- 
up which will lead to production of pure SC. The spin con- 
ductance is defined as G'?(G'?) oc \t\^ + I^aP whereas the 



given by 



In(fcd) where a = (52 — 29i)/27rfi,Wi^, d is 



charge conductance is given by G'^{G'^) cx — (|tp — I^Ap)- 
The t and | arrows in the subscript represent the spin po- 
larization of the injected electrons from the ferromagnetic 
lead (see Fig. 1). The negative sign in the expression for 
G'i^(G^) arises because it is a sum of contribution com- 
ing from two oppositely charged particles (electrons and 
holes). The first S'-matrix, has r = (reflection- less) , 
r A 7^ and t = tA- This is not a fixed point and hence the 
parameters of the S'-matrix will flow under RG . It is easy 
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Fig. 3; The variation of — — I^Ap) <x. G'^ or and the variation of + |tAp) G| or G| are plotted as a function of the 
dimcnsionless parameter / where I = ln(L/d) and L is either = hvp/kgT at zero bias or Ly = hvp/eV at zero temperature and d is 
the short distance cut-off for the RG fiow. The three curves in each plot correspond to three different values of V{0) and V(2kp) for the 
NSN junction. These plots correspond to the S-matrix given by Sa- 



to see from Eqs. 4-7, that for this case, the RG equations 
for t and tj^ are identical, and hence it is ensured that the 
RG flow will retain the equality of the t and leading to 
the preservation of pure SC . Physically this implies that 
if we start the experiment with this given S'-matrix (Si) 
at the high energy scale (at finite bias voltage and zero 
temperature or at zero bias and finite temperature), then, 
as we reduce the bias in the zero temperature case (or 
reduce the temperature in the zero bias case), the correla- 
tions arising due to inter-electron interactions in the wire 
are such that the amplitude of t and tA will remain equal 
to each other. The quantity which increases with increas- 
ing length scale L is the absolute value of the amplitude 
t or tA leading to a monotonic increase of pure SC till it 
saturates at the maximum value allowed by the symme- 
tries of the S'-matrix, Si (Fig. 2). Here all the S'-matrix 
elements are assumed to be energy independent and hence 
the bias dependence is solely due to RG flow. Of course 
the bias window has to be small enough so that the en- 
ergy dependence of t, tA, r and ta can be safely ignored. 
This saturation point is actually a stable fixed point of 
the theory if the junction remains reflection-less. So we 
observe that the transmission (both t and t^) increases 
to maximum value while the AR amplitude scales down 
to zero. This flow direction is quite different from that 
of the standard case of a single impurity in an interacting 
electron gas in 1~D where any small but finite reflection 
amplitude gets enhanced under RG flow ultimately leading 
to zero transmission [16]. The difference here is because 
the RG flow is solely due to the existence of the finite pair 



potential (due to r^) and not due to the usual Friedel os- 
cillations (due to r). Hence the electrons in the wire have 
an effective attractive interaction leading to a counter in- 
tuitive RG flow. We remark that the interaction induced 
correction enhances the amplitude for pure SC and also 
stabilizes the pure SC operating point. This makes the 
operating point, Si quite well-suited for an experimental 
situation. Fig. 2 shows the variation of the pure spin con- 
ductance (= 2 X |tp in units of e'^/h) as a function of 
relevant length scale, L of the problem. 

The second case corresponds to the most symmetric S- 
matrix (S'2). It is a fixed point of RG equations and is given 
by r = 1/2, TA = ~l/2,t = 1/2, tA = 1/2. Here also t is 
equal to tA as in the previous case and thus the junction 
will act like a perfect charge filter resulting in pure SC in 
the right wire (if spin polarized charge current is injected 
in the left wire). However, this S'-matrix (S'2) represents 
an unstable fixed point. Due to any small perturbation, 
the parameters tend to fiow away from this unstable fixed 
point to the most stable disconnected fixed point given by 
\r\ = 1 as a result of RG fiow. So this S'-matrix (S'2) is 
not a stable operating point for the production of pure SC. 
But, it is interesting to note that if we switch on a small 
perturbation around this fixed point, the charge conduc- 
tance exhibits a non-monotonic behavior under RG fiow 
(Fig. 3). This non-monotonicity results from two compet- 
ing effects viz., transport through both electron and hole 
channels and, the RG fiow of gi, g2- This essentially leads 
to negative differential conductance (N DC) [17]. Elaborat- 
ing it further, all it means is that if we start an experiment 
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Fig. 4; The variation of — — ItAp) oc G'^ or Gj^ and the variation of + |tAp) or are plotted in left and right panel 

plots as a function of the dimensionless parameter I where I = ln{L/d) and L is either Lx = hvp/kgT at zero bias or Ly = hvp/eV at 
zero temperature and d is the short distance cut-off for the RG flow. The three curves in each plot correspond to three different values of 
y(0) and V{2kp) for the FSN junction. These plots correspond to the S-matrix given by S3. 



with this given S'-matrix (S2) at zero temperature and at 
finite bias, then as we go towards zero bias, the conduc- 
tance will show a rise with decreasing bias for a certain 
bias window. This can be seen from Fig. 3. This aspect 
of the RG flow can be of direct relevance for manipulating 
electron and spin transport in some mesoscopic devices. 

Now we will switch to the case of ferromagnetic 
half metal— superconductor— normal metal (FSN) junc- 
tion which comprises of a 1-D ferromagnetic half metal 
(assuming | polarization) on one side and a nor- 
mal 1-D metal on the other side (in a way similar 
to the set-up shown in Fig. 1). This case is very 
complicated to study theoretically because the mini- 
mal number of independent complex-valued parameters 
that are required to parameterize the S'-matrix is nine 
as opposed to the previous (symmetric) case which 
had only four such parameters. These are given by 

„11 „22 22 ^12 J.2I 22 22 ^12 J ^21 TT 

'•fp '-f-f) 'lit '■AfT' All' -4tT' All' TT' O'^-''-' "-ff- -Tieie, 
1(2) is the wire index for the ferromagnetic (normal) wire 
while, t and | are the respective spin polarization in- 
dices for the electron. The large number of independent 
parameters in this case arise because of the presence of 
ferromagnetic half-metallic wire which destroys both the 
spin rotation symmetry and the left-right symmetry. The 
only remaining symmetry is the particle-hole symmetry. 
Analogous to the RG equations (given by Eqs. 4-7) for 
the NSN case, it is possible to write down all the nine 
RG equations for FSN case and solve them numerically to 
obtain the results as shown in Fig. 4. (For further details, 
see Ref. [18]). In this case, the elements of a represen- 



tative S'-matrix (S3) which correspond to the production 
of pure SC are |r|}| = \rf^\ = jr^J] = \t\^^^\ = \t\\^\ = 

kin' = l^'^ul = I^TtI = I^TtI = I/a/S and the correspond- 
ing phases associated with each of these amplitudes are 
7r/3, TT, 0, — 7r/3, 0, 7r/3, 0, tt, — 7r/3 respectively. By solving 
the nine coupled RG equations for the above mentioned 
nine independent parameters, we have checked numeri- 
cally that this is not a fixed point of the RG equation and 
hence it will flow under RG and finally reach the trivial 
stable fixed point given by r|| = r^^ = rJJ = 1. Now if we 
impose a bias on the system from left to right, it will create 
a pure SC on the right wire because \t]l^^ \ is exactly equal 
to But of course, this is a highly unstable operating 

point for production of pure SC as this is not even a fixed 
point and hence will always flow under any variation of 
temperature or bias destroying the production of pure SC. 
In this case also, the spin conductance shows a monotonic 
behavior while, the charge conductance is non-monotonic 
and hence will have NDC in some parameter regime. It 
is worth noticing that in this case the interaction param- 
eters gi and 32 both do not scale on the left wire as it is 
completely spin polarized while gi and 32 do scale on the 
right wire as it is not spin polarized. Hence even if we be- 
gin our RG flow with symmetric interaction strengths on 
both left and right wires, they will develop an asymmetry 
under the RG flow. 

Finally, we consider another important aspect 
that nicely characterizes these hybrid structures 
from a spintronics application point of view. If 
the QW on the two sides of the superconductor are 
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ferromagnetic half metals then we have a junction 
of ferromagnet— superconductor— ferromagnet (FSF). 
We calculate the tunnelling magnctorcsistance ratio 
(TMR ) [2] which is defined as follows 



TMR 



G|| — G 



n 



G 



n 



(8) 



Here, G|| corresponds to the conductance across the junc- 
tion when both left and right wires are in parallel spin- 
polarized configurations. G-\i corresponds to the case 
when the left and right wires are in anti-parallel spin- 
polarized configurations. Thus, TMR is the maximum rel- 
ative change in resistance in going from the parallel to 
the anti-parallel configuration. For the parallel case, the 
CAR amplitude [tA) is zero and the only process which 
contributes to the conductance is the direct tunnelling 
process. This is because the CAR process involves non- 
local pairing of "f e in the left wire with j e in the right 
wire to form a Cooper pair. However for [ e , the den- 
sity of states is zero in the right wire which makes this 
process completely forbidden. Hence, G|| ex \t\^ . On 
the other hand, for the anti-parallel case, G|| oc — I^aP 
as there is no density of states for the t e in the right 
lead and so no direct tunnelling of t e across the junc- 
tion is allowed; hence CAR is the only allowed process. 
Note that the negative sign in G|| leads to a very large 
enhancement of TMR (as opposed to the case of standard 
ferromagnet— normal metal— ferromagnet (FNF) junction) 
since the two contributions will add up. A related set-up 
has been studied in [19] where also a large TMR has been 
obtained. 

One can then do the RG analysis for both parallel and 
anti-parallel cases. It turns out that the equations for 
\t\ and \tA\ are identical leading to identical temperature 
(bias) dependance. The RG equation for \tA\ is 



1f = -^'- [i-M' 



(9) 



Here, /? = {g2 — gi)/2-KhvF- \t\ satisfies the same equa- 
tion. So, in a situation where the reflection amplitudes at 
the junction for the two cases are taken to be equal then 
it follows from Eq. 8 that the TMR will be pinned to its 
maximum value i.e. magnitude of TMR ~ 2 and the tem- 
perature dependence will be flat even in the presence of 
intcr-elcctron interactions. 

Conclusions. — In this letter, we have studied both 
spin and charge transport in NSN, FSN, and FSF struc- 
tures in the context of 1-D QW. We calculated the cor- 
rections to spin and charge transport arising from inter- 
electron interactions in the QW. We demonstrated the pos- 
sibility for production of pure SC in such hybrid junctions 
and analysed its stability against temperature and volt- 
age variations. Finally, we also showed that the presence 
of the CAR process heavily enhances the TMR in such ge- 
ometries. 
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Abstract. - We demonstrate possible scenarios for production of pure spin current and large tun- 
nelling magnetoresistance ratios from elastic co-tunnelling and crossed Andreev reflection across 
a superconducting junction comprising of normal metal-superconductor-normal metal, where, the 
normal metal is a one-dimensional interacting quantum wire. We show that there are fixed points 
in the theory which correspond to the case of pure spin current. We analyze the influence of 
electron-electron interaction and see how it stabilizes or de-stabilizes the production of pure spin 
current. These fixed points can be of direct experimental relevance for spintronics application 
of normal metal-superconductor-normal metal junctions of one- dimensional quantum wires. We 
also calculate the power law temperature dependence of the crossed Andreev reflection enhanced 
tunnelling magnetoresistance ratio for the normal metal-superconductor-normal metal junction. 
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. Introduction. — Two fundamental degrees of free- 
dom associated with an electron that are of direct interest 
to condensed matter physics are its charge and spin. Until 
V^ery recently, all conventional clcctron-bascd devices have 
been solely based upon the utilization and manipulation 
of the charge degree of freedom of an electron. However, 
the realization of the fact that devices based on the spin 
degree of freedom can be almost dissipation-less and with 
very fast switching times, has led to an upsurge in re- 
search activity in this direction in recent years [?,?,?]. 
The first step towards realization of spin-based electronics 
(spintronics) would be to produce pure spin current (SC). 
From a purely theoretical point of view, it is straight for- 
ward to define a charge current as a product of local charge 
density with the charge velocity, but such a definition can- 
not be straight forwardly extended to the case of SC. This 
is because both spin S and velocity v are vector quantities 
and hence the product of two such vectors will be a tensor. 

In this letter, we adopt the simple minded definition 
of SC, which is commonly used [?]. It is just the prod- 
uct of the local spin polarization density associated with 
the electron or hole, (a scalar s which is either positive 
for up-spin or negative for down-spin) and its velocity [?]. 
The most obvious scenario in which one can generate a 
pure SC in the sense defined above would be to have (a) 



an equal and opposite flow of identically spin-polarized 
electrons through a channel, such that the net charge 
current through the channel is nullified leaving behind a 
pure SC, or (b) alternatively, an equal fiow of identically 
spin polarized electrons and holes in the same direction 
through a channel giving rise to pure SC with perfect can- 
cellation of charge current. In this letter, we explore the 
second possibility for generating pure SC using a normal 
metal— superconductor— normal metal (NSN) junction. 

Proposed device and its theoretical modelling. — 

The configuration we have in mind for the production of 
pure SC is shown in Fig. 1. The idea is to induce a pair 
potential in a small region of a quantum wire (QW) by de- 
positing a superconducting strip on top of the wire (which 
may be, for instance, a carbon nanotube) due to proximity 
effects. If the strip width on the wire is of the order of the 
phase coherence length of the superconductor, then both 
direct electron to electron co-tunnelling as well as crossed 
Andreev electron to hole tunnelling can occur across the 
superconducting region [?]. It is worth pointing out that 
in the case of a singlet superconductor, which is the case 
we consider, both the tunnelling processes will conserve 
spin. In order to describe the mode of operation of the 
device (see Fig. 1), we first assume that the S'-matrix rep- 
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Fig. 1: A 1-D quantum wire (carbon nanotubc) connected to a ferromagnetic (F) lead on the left and a normal (N) lead on the right. 
The thin strip in the middle of the wire depicts a 2-D layer of superconducting material deposited on top of the wire. 



resenting the NSN junction described above respects par- 
ity symmetry about the junction, particle-hole symmetry 
and spin-rotation symmetry. Considering all the symme- 
tries, we can describe the superconducting junction con- 
necting the two half wires by an S'-matrix with only four 
independent parameters namely, (i) the normal refection 
amplitude (r) for e (h), (ii) the transmission amplitude [t) 
for e (h), (Hi) the Andreev reflection (AR) amplitude (r^) 
for e (h), and (iv) the crossed Andreev reflection (CAR) 
amplitude {Ia) for e (h). If we inject spin polarized elec- 
tron (t e) from the left QW using a ferromagnetic contact 
and tune the junction parameters such that t and ta are 
equal to each other, it will lead to a pure SC flowing in 
the right QW (see Fig. 1). This is so because, on an av- 
erage, an equal number of electrons (t e) (direct electron 
to electron tunnelling) and holes (t h) (crossed Andreev 
electron to hole tunnelling) are injected from the left wire 
to the right wire resulting in production of pure SC in the 
right wire. Note that spin up holes (t h) implies a Fermi 
sea with an absence of spin down electron (which is what 
is needed for the incident electron (t e) to form a Cooper 
pair and enter the singlet superconductor). 

Renormalisation Group study. Wc now include 
the effects of inter-electron interactions on the S'-matrix 
using the renormalisation group (RG) method introduced 
in Ref . [?] . This was furthur generalized to the case of mul- 
tiple wires in Refs. [?,?] and to the case of 1-D NS junc- 
tion [?, ?, ?]. The basic idea of the method is as fol- 
lows. The presence of scattering (reflection) induces 
Friedel oscillations in the density of non-interacting elec- 
trons. Within a mean field picture for a weakly inter- 
acting electron gas, the electron not only scatters off the 
potential barrier but also scatters off these density oscil- 
lations with an amplitude proportional to the interaction 
strength. Hence by calculating the total reflection ampli- 
tude due to scattering from the scalar potential scatterer 
and from the Friedel oscillations created by the scatterer, 
we can include the effect of electron-electron interaction in 
calculating transport. This approach can be generalized 



to junctions of one-dimensional (1-D) QW with a single 
superconductor. In this case, there will be non-zero AR in 
the bulk of the wire due to proximity induced pair poten- 
tial, besides the AR right at the NSN junction which turns 
an incoming electron into an outgoing hole. 

The RG equations for a NS case have been obtained 
earlier using bosonization [?,?,?] and using WIRG [?,?]. In 
this letter wc extend these WIRG results to the NSN case. 
To obtain the RG equations in the presence of AR and 
CAR for the NSN junction, we follow a procedure similar to 
that used in Ref. [?] . The fcrmion fields expanded around 
the left and right Fermi momenta on each wire can be 
written as, V»s(x) = (a;) '^^ " ~f ■^o^s{x)e-'^^'' ; 
where i is the wire index, s is the spin index which can be 
t,i and I{0) stand for incoming (outgoing) fields. Note 
that ^'/(o)(x) arc slowly varying fields on the scale of fc^^^. 
For a momentum in the vicinity of fc^?, the incoming and 
outgoing fields can be Fourier expanded as: 



*fcs(a;) 



dk 



gj (fe + ^f) a; _|_ ^i{-k + kp) X 



(1) 



where 6^,^ is the electron destruction operator and d^^ is 
the hole destruction operator and we have allowed for non- 
conservation of charge due to the proximity effect. We 
allow for short-range density-density interactions between 

the fermions, Hint ~ 2 I dx dy p^gV{x — y) p^^, where the 
sum over the spin indices is assumed. 

Then the effective Hamiltonian, can be derived using a 
Hartree— Fock (HF) decomposition of the interaction. The 
charge conserving HF decomposition leads to the interac- 
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Fig. 2: The variation of |tp {=|iAp) is plotted as a function of 
the dimensionless parameter / where I = ln(L/d) and L is either 
Lx = hvp/kgT at zero bias or Ly = hvp/eV at zero temperature 
and d is the short distance cut-off for the RG fiow. The three curves 
correspond to three different values of V{0) and V{2kp) for the 
NSN junction. This case corresponds to the 5-matrix given by Si. 



tion Hamiltonian (normal) of the form 
-i{g2 - 2gi) f°° dx 



njN 

'lint — 



47r 



a 



(2) 



(We have assumed spin-symmetry and used r^i = = 
Ti.) This has been derived earlier in Ref. [?]. We use the 
same method, but now, we also allow for a charge non- 
conserving HF decomposition and arrive at the (Andreev) 
Hamiltonian 



ri\nt — 



-*(.9i + .92) 



dx 



(3) 



Note that, even if this appears non-charge conserving, 
charge conservation is taken care of by the 2e charge that 
flows into the superconductor due to proximity induced 
Cooper pairing. 

In Ref. [?], the perturbatively calculated correction to 
the reflection amplitude under exp [—iTi.^^ t] was derived 
to first order in a. For electrons with spin incident with 
momentum k with respect to kp, this was shown to be 



given by 



In(fcd) where a = (52 — 29i)/27rfi,Wi^, d is 



a short distance cutoff and gi = V{2kp), g2 = V{0). It 
is worth noting that both gi and (72 scale under RG [?]. 
Hence, the values given above for gi and 172 are the micro- 
scopic values it takes in the original Hamiltonian defined 
at the short distance scale. 

Analogously, here, we calculate the amplitude to go 
from an incoming electron wave to an outgoing hole wave 

a' TAi 



under exp [— zTij^^ t]. It is given by 



■ In(fcd) where 



Q^' ~ {91 + g2)/2TThvp. Note that the spin of the outgo- 
ing hole is always the same as the spin of the incoming 
electron, since the Andreev Hamiltonian also conserves 
spin for a singlet superconductor. We see that there is a 
logarithmic singularity at the fc ^ limit which implies 
that the lowest order perturbation theory is not enough to 
calculate correction to the reflection and AR amplitudes 
when the momenta of the incident particles are very close 
to the Fermi wave vector. Following Yue et al. [?], we sum 
up these most divergent processes using the "poor man's 
scaling" approach [?] to obtain RG equations for the nor- 
mal refection amplitude (r), the transmission amplitude 
(i), the AR amplitude (r^), and the CAR amplitude (i^) 
which are as follows, 



dr 

m 

drA 
dl 

dt 
Jl 

dtA 
dl 



r\+t\) -r{l-\r\^)) 



~a' {r|rAp + r\tAt] 

- a {\r\^rA + ttAr*] + y {ta 

- a {\r\^t + r*rAtA] 
-a' {\rA\h + rr\tA] 

- a {rV^i + \r\'^tA] 
-a {rtr\ + \rA\^tA] 



(4) 



(5) 



(6) 



(7) 



Note that when < = = these equations reduce to the 
RG equations obtained in Ref. [?] for the case of NS junc- 
tion. 

Results and Discussion. — We propose two possible 
S*- matrices {Si and S2) that can be realized within our set- 
up which will lead to production of pure SC. The spin con- 
ductance is defined as G^{G^^) oc \t\^ + I^aP whereas the 
charge conductance is given by Gf'(Gf) cx —{\t\^ — I^Ap)- 



The t and | arrows in the subscript represent the spin po- 
larization of the injected electrons from the ferromagnetic 
lead (see Fig. 1). The negative sign in the expression for 
G|^(G^) arises because it is a sum of contribution com- 
ing from two oppositely charged particles (electrons and 
holes). The first S'-matrix, has r = (reflection- less) , 
r A 7^ and t = tA- This is not a fixed point and hence the 
parameters of the S'-matrix will flow under RG . It is easy 
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Fig. 3; The variation of — — I^Ap) <x. G'^ or and the variation of + |tAp) G| or G| are plotted as a function of the 
dimcnsionless parameter / where I = ln(L/d) and L is either = hvp/kgT at zero bias or Ly = hvp/eV at zero temperature and d is 
the short distance cut-off for the RG fiow. The three curves in each plot correspond to three different values of V{0) and V(2kp) for the 
NSN junction. These plots correspond to the S-matrix given by Sa- 



to see from Eqs. 4-7, that for this case, the RG equations 
for t and tj^ are identical, and hence it is ensured that the 
RG flow will retain the equality of the t and leading to 
the preservation of pure SC . Physically this implies that 
if we start the experiment with this given S'-matrix (Si) 
at the high energy scale (at finite bias voltage and zero 
temperature or at zero bias and finite temperature), then, 
as we reduce the bias in the zero temperature case (or 
reduce the temperature in the zero bias case), the correla- 
tions arising due to inter-electron interactions in the wire 
are such that the amplitude of t and tA will remain equal 
to each other. The quantity which increases with increas- 
ing length scale L is the absolute value of the amplitude 
t or tA leading to a monotonic increase of pure SC till it 
saturates at the maximum value allowed by the symme- 
tries of the S'-matrix, Si (Fig. 2). Here all the S'-matrix 
elements are assumed to be energy independent and hence 
the bias dependence is solely due to RG flow. Of course 
the bias window has to be small enough so that the en- 
ergy dependence of t, tA, r and ta can be safely ignored. 
This saturation point is actually a stable fixed point of 
the theory if the junction remains reflection-less. So we 
observe that the transmission (both t and t^) increases 
to maximum value while the AR amplitude scales down 
to zero. This flow direction is quite different from that 
of the standard case of a single impurity in an interacting 
electron gas in 1-D where any small but finite reflection 
amplitude gets enhanced under RG flow ultimately lead- 
ing to zero transmission [?] . The difference here is because 
the RG flow is solely due to the existence of the finite pair 



potential (due to r^) and not due to the usual Friedel os- 
cillations (due to r). Hence the electrons in the wire have 
an effective attractive interaction leading to a counter in- 
tuitive RG flow. We remark that the interaction induced 
correction enhances the amplitude for pure SC and also 
stabilizes the pure SC operating point. This makes the 
operating point, Si quite well-suited for an experimental 
situation. Fig. 2 shows the variation of the pure spin con- 
ductance (= 2 X |tp in units of e'^/h) as a function of 
relevant length scale, L of the problem. 

The second case corresponds to the most symmetric S- 
matrix (S2). It is a fixed point of RG equations and is given 
by r = 1/2, = -1/2, t = 1/2, tA = 1/2. Here also t is 
equal to tA as in the previous case and thus the junction 
will act like a perfect charge filter resulting in pure SC in 
the right wire (if spin polarized charge current is injected 
in the left wire). However, this S'-matrix (S'2) represents 
an unstable fixed point. Due to any small perturbation, 
the parameters tend to fiow away from this unstable fixed 
point to the most stable disconnected fixed point given by 
\r\ = 1 as a result of RG fiow. So this S'-matrix (S'2) is 
not a stable operating point for the production of pure SC. 
But, it is interesting to note that if we switch on a small 
perturbation around this fixed point, the charge conduc- 
tance exhibits a non-monotonic behavior under RG fiow 
(Fig. 3). This non-monotonicity results from two compet- 
ing effects viz., transport through both electron and hole 
channels and, the RG flow of gi, 172. This essentially leads 
to negative differential conductance (NDC) [?]. Elaborat- 
ing it further, all it means is that if we start an experiment 
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Fig. 4; The variation of — — ItAp) oc G'^ or Gj^ and the variation of + |tAp) or are plotted in left and right panel 

plots as a function of the dimensionless parameter I where I = ln{L/d) and L is either Lx = hvp/kgT at zero bias or Ly = hvp/eV at 
zero temperature and d is the short distance cut-off for the RG flow. The three curves in each plot correspond to three different values of 
y(0) and V{2kp) for the FSN junction. These plots correspond to the S-matrix given by S3. 



with this given S'-matrix (S2) at zero temperature and at 
finite bias, then as we go towards zero bias, the conduc- 
tance will show a rise with decreasing bias for a certain 
bias window. This can be seen from Fig. 3. This aspect 
of the RG flow can be of direct relevance for manipulating 
electron and spin transport in some mesoscopic devices. 

Now we will switch to the case of ferromagnetic 
half metal— superconductor— normal metal (FSN) junc- 
tion which comprises of a 1-D ferromagnetic half metal 
(assuming | polarization) on one side and a nor- 
mal 1-D metal on the other side (in a way similar 
to the set-up shown in Fig. 1). This case is very 
complicated to study theoretically because the mini- 
mal number of independent complex-valued parameters 
that are required to parameterize the S'-matrix is nine 
as opposed to the previous (symmetric) case which 
had only four such parameters. These are given by 

„11 „22 22 ^12 J.2I 22 22 ^12 J ^21 TT 

'•fp '-f-f) 'lit '■AfT' All' -4tT' All' TT' O'^-''-' "-ff- -Tieie, 
1(2) is the wire index for the ferromagnetic (normal) wire 
while, t and | are the respective spin polarization in- 
dices for the electron. The large number of independent 
parameters in this case arise because of the presence of 
ferromagnetic half-metallic wire which destroys both the 
spin rotation symmetry and the left-right symmetry. The 
only remaining symmetry is the particle-hole symmetry. 
Analogous to the RG equations (given by Eqs. 4-7) for 
the NSN case, it is possible to write down all the nine 
RG equations for FSN case and solve them numerically to 
obtain the results as shown in Fig. 4. (For further details, 
see Ref. [?]). In this case, the elements of a represen- 



tative S'-matrix (S3) which correspond to the production 
of pure SC are |r|}| = \rf^\ = jr^J] = \t\^^^\ = \t\\^\ = 

kin' = l^'^ul = I^TtI = I^TtI = I/a/S and the correspond- 
ing phases associated with each of these amplitudes are 
7r/3, TT, 0, — 7r/3, 0, 7r/3, 0, tt, — 7r/3 respectively. By solving 
the nine coupled RG equations for the above mentioned 
nine independent parameters, we have checked numeri- 
cally that this is not a fixed point of the RG equation and 
hence it will flow under RG and finally reach the trivial 
stable fixed point given by r|| = r^^ = rJJ = 1. Now if we 
impose a bias on the system from left to right, it will create 
a pure SC on the right wire because \t]l^^ \ is exactly equal 
to But of course, this is a highly unstable operating 

point for production of pure SC as this is not even a fixed 
point and hence will always flow under any variation of 
temperature or bias destroying the production of pure SC. 
In this case also, the spin conductance shows a monotonic 
behavior while, the charge conductance is non-monotonic 
and hence will have NDC in some parameter regime. It 
is worth noticing that in this case the interaction param- 
eters gi and 32 both do not scale on the left wire as it is 
completely spin polarized while gi and 32 do scale on the 
right wire as it is not spin polarized. Hence even if we be- 
gin our RG flow with symmetric interaction strengths on 
both left and right wires, they will develop an asymmetry 
under the RG flow. 

Finally, we consider another important aspect 
that nicely characterizes these hybrid structures 
from a spintronics application point of view. If 
the QW on the two sides of the superconductor are 
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ferromagnetic half metals then we have a junction We acknowledge use of the Beowulf cluster at H.R.I, 
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We calculate the tunnelling magnctorcsistance ratio gramme at WIS, Israel. 
(TMR ) [?] which is defined as follows 



TMR 



G^^ - G^i 



G 



n 



(8) 



Here, G|| corresponds to the conductance across the junc- 
tion when both left and right wires are in parallel spin- 
polarized configurations. G-\i corresponds to the case 
when the left and right wires are in anti-parallel spin- 
polarized configurations. Thus, TMR is the maximum rel- 
ative change in resistance in going from the parallel to 
the anti-parallel configuration. For the parallel case, the 
CAR amplitude [tA) is zero and the only process which 
contributes to the conductance is the direct tunnelling 
process. This is because the CAR process involves non- 
local pairing of "f e in the left wire with j e in the right 
wire to form a Cooper pair. However for [ e , the den- 
sity of states is zero in the right wire which makes this 
process completely forbidden. Hence, ex \t\^ . On 
the other hand, for the anti-parallel case, G|| (x — I^aP 
as there is no density of states for the t e in the right 
lead and so no direct tunnelling of t e across the junc- 
tion is allowed; hence CAR is the only allowed process. 
Note that the negative sign in G|| leads to a very large 
enhancement of TMR (as opposed to the case of standard 
ferromagnet— normal metal— ferromagnet (FNF) junction) 
since the two contributions will add up. A related set-up 
has been studied in [?] where also a large TMR has been 
obtained. 

One can then do the RG analysis for both parallel and 
anti-parallel cases. It turns out that the equations for 
\t\ and \tA\ are identical leading to identical temperature 
(bias) dependance. The RG equation for \tA\ is 

^^ = -PtA [l-\tA?] (9) 

Here, /? = (,92 — gij/'^T^^vp- \t\ satisfies the same equa- 
tion. So, in a situation where the reflection amplitudes at 
the junction for the two cases are taken to be equal then 
it follows from Eq. 8 that the TMR will be pinned to its 
maximum value i.e. magnitude of TMR —2 and the tem- 
perature dependence will be flat even in the presence of 
intcr-elcctron interactions. 

Conclusions. — In this letter, we have studied both 
spin and charge transport in NSN, FSN, and FSF struc- 
tures in the context of 1-D QW. We calculated the cor- 
rections to spin and charge transport arising from inter- 
electron interactions in the QW. We demonstrated the pos- 
sibility for production of pure SC in such hybrid junctions 
and analysed its stability against temperature and volt- 
age variations. Finally, we also showed that the presence 
of the CAR process heavily enhances the TMR in such ge- 
ometries. 



